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(57) Abstract 

Apparatus for generating a diagnostic display of a patient's seismocardiogram (SCG). The SCG waveform is collected (50) 
and segmented into wavelets (56). The wavelets are cross-correlated (60-65) and averaged to generate a representative waveform 
(64). Certain measurements are performed on the waveform and these parameters (68) are multiplied by selected weighting coeffi- 
cients to generate an advisory diagnosis (76). 
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SEISMOCARDIOGRAPHIC ANALYSIS SYSTEM 



Background of the Invention 
Field of the Invention 
5 The present invention relates to seisznocardiography , 

and more particularly to a non- invasive seismocar- 
diograph system for diagnosing abnormal cardiac func- 
tion. The system collects, categorizes, averages and 
analyses the patient's seismocardiogram (SCG) and pro- 
10 vides the attending physician with a preliminary and 

advisory diagnosis based upon a secpience of measurements 
made upon the collected data. 
Description of the Background Art 

Seismocaxdiography is a technique for the non- 
15 invasive assessment of the mechanical motion of the 

heart. In general, compress ional waves generated by the 
mechanical motion of the heart are detected by a seismic 
transducer from the surface of the patient's chest. The 
preferred location for the transducer is typically the 
20 lower third of the sternum. 

U.S. Patent No. 4,989,611 to Zanetti, et al . , issued 
Feb. 5, 1991 discloses a system for the non-invasive 
assessment of myocardial ischemia through the use of 
seismocardiography . This patent teaches, inter alia, 
25 the exclusion of electrically abnoirmal beats from the 
collection of SCG data used for diagnostic purposes . 

Although seismocardiography has proven to be a 
reliable and accurate method of assessing the mechanical 
function of the heart, it has not achieved dominance as 
30 a non-invasive diagnostic procedure, in part due to the 
unf amiliarity with appropriate interpretive technique 
among physicians /diagnosticians, 

Suimnarv of the Invention 
In contrast to prior seismocardiographic systems , 
35 the present invention formulates a preliminary diagnosis 
(78) based upon certain measurements made upon collected 
and processed data. 

In operation, patient data is taken under a stress 
test protocol (50). Simultaneous electrocardiographic 
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ECG (20) and seismocardiographic SCG (10) data are taken 
pre-exerciser post-exercise and dtiring recovery from 
exercise . 

The SCG data is broken up into discrete waveform 
5 segments (56) which correspond approximately to the 

mechanical motion of the heart during individual heart- 
beats. Each of these SCG waveform segments is called a 
"wavelet " ( wl-w6 ) . 

Each SCG wavelet is compared with all other wavelets 

10 to ascertain their degree of similarity (58). This 
cross-correlation process is repeated iteratively. 

The cross -correlation process categorizes and groups 
wavelets into families (a-e) based upon intra- family 
similarity. It has been found through experimentation 

15 that as many as ninety percent (90%) of the wavelets 

will fall into one family (c) . A family of highly self- 
similar wavelets^ which contains the majority of all 
wavelets analyzed is formed and is referred to as the 
"dominant family of wavelets" (64). 

20 Families of dominant wavelets are extracted for two 

phases of the stress test. A pre-exercise or baseline 
dominant family is developed as is a post-exercise 
dominant family. 

Each dominant family set of wavelets is averaged to 

25 arrive at a canonical form wavelet which typifies the 
data collected during the corresponding phase of the 
stress test. This canonical form wavelet is referred to 
as the "representative waveform". The representative 
waveforms generated are displayed to the diagnostician 

30 (see 82 and 83rFIG.3 and 84 on FIG. 5). Two of the three 
phases of the stress test are used to generate the 
advisoiY diagnosis . 

Certain features or points on the representative 
waveforms can be consistently identified (66). These 

35 point are identified on the waveforms and certain 
measurements are. made between points . For example 
amplitudes are measured and expressed in acceleration 
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unlt:s (g). Slopes are measured and are expressed as 
accelera-tlon per unit time ( g\ms ) . Time interval between 
events are also made and expressed in milliseconds (ms) • 

This measurement and calculation process is referred 
5 to as "parameterization" of the representative waveform 
(68). The resultant times, amplitudes, slopes and ratios 
are called collectively "parameters", and these are 
displayed in a hard copy format as shown in FIG. 4. 

These parameters are multiplied by certain pre- 
10 defined "weighting coefficients" to form "weighted para- 
meters" (70) . The weighted parameters are added up and 
if their collected value tends toward zero, the can- 
didate waveform is likely to have arisen from a healthy 
heart, while if the sum tends toward unity then the 
15 candidate waveform is indicative of a disease state 
(74). 

The advisory diagnosis thus formed, is presented 78 
to the diagnostician along with other waveform and para- 
meter information to aid in interpretation by the 
20 diagnostician. 

Brief Description of the Drawings 
Throughout the figures, identical reference 
numerals refer to identical structural elements, 
wherein: 

25 FIG. 1 is a schematic view depicting the 

seismocardiographic analysis system; 

FIG. 2 is a schematic plan view of the dia- 
gnostic waveforms indicating SCG nomenclature; 

FIG. 3 is a panel of data in the foannat 
30 presented to the diagnostician; 

FIG. 4 is a panel of data in the format pre- 
sented to the diagnostician; 

FIG. 5 is a panel of data in the format pre- 
sented to the diagnostician; 
35 FIG. 6 is a schematic presentation showing the 

segmentation of the SCG data into wavelets; 
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FIG. 7 is a histogram showing the results of a 
hypothetical cross-correlation of the schematic data of 
FIG. 6; 

PIG. 8 is a table of parameters along with 
5 their identifier code and units of measure; 

FIG. 9 is a tedale of parameter weighting 
coefficients and their identifier code, as computed by a 
neural network, presenting values for pre-exercise or 
baseline; post-exercise, and recovery-from-exercise 
10 phases of the stress test value, and; 

FIG. 10 is a flow chart depicting program flow 
for the seismocardiographic analysis system of FIG. 1. 

Detailed Description of the Invention 
Overview of the Physical System 
15 A wide range of hardware platforms can be used to 

carry out the invention. FIG. 1 shows an exemplary 
system which includes a seismocardiographic sensor 10, 
placed on the chest of the patient 12 to collect SCG 
data from the patient. In a preferred mode, electrocar- 
20 diographic data, from ECG leads 20, is collected and 
supplied to the seismocardiograph 14. 

The seismocardiograph 14 itself, is partitioned into 
computing apparatus 22, signal processing hardware 16, 
and display structures 18. The computing apparatus 
25 collects and processes the SCG and ECG data while the 
display 18 provide graphic output to the physician in 
both hard copy and display screen formats. FIG. 3, 4 and 5 
represent hard copy output provided by the machine to 
the diagnostician. 
30 The preferred seismocardiographic sensor is a ultra- 

low frequency seismic accelerometer which weighs about 
one kilogram and which may be conveniently placed on the 
patients sternum. 

Any one of a number of suitable ECG electrode arrays 
35 may be used to carry out the invention. 

Both the SCG and ECG signals are buffered and 
isolated and may be filtered using a bandpass filter 
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having a .3 -to 50 Hz. passband. The frequency response 
should be ex-tended t:oifard DC as much as possible consis- 
tent with coupling requirements. Although the high fre- 
quency cut off Is relatively low, much higher frequency 
5 response for the sensor Is desired to preserve the 

fidelity of the SC6 waveform. The SCG and EC6 data Is 
then converted to digital form at a sampling rate of 250 
sample /second. An IBU compatible AT class computer 22 
Is suitable for data reduction and processing. 
10 Overview of System Operation 

In general, the analog ECG and SCG signals are 
collected from the patient through the signal processing 
module 16 . This data is converted into a digital format 
and supplied to the computing module 22 . This activity 
15 corresponds to process 50 shown on FIG. 10. 

The first objective of the computing system 16 is to 
parse the SCG data based on the ECG data and segment the 
SCG waveform data into wavelets. This object requires 
elimination of non-sinus beats and the identification of 
20 the R-wave peak in the remaining waveform data. This 
process corresponds to process 52 on FIG. 10. 

As shown in FIG. 6 the ECG waveform 24 is scanned 
and an appropriate collection of rules is applied to 
determine the locations of R-wave peaks 1 through 6 in 
25 the ECG channel data. 

It is Important to exclude SCG data arising from 
non-sinus heartbeats. It has been determined that non- 
sinus beats are so abnormal mechanically, that they 
degrade the diagnostic potential of the SCG and should 
30 be considered noise. Therefore the ECG channel is 

scanned to exclude all non-sinus beats and to determine 
the location of the remaining R waves. A variety of ECG 
beat classification schemes may be used without 
departing from the scope and spirit of the invention. It 
35 is important that the ECG classifier algorithm be con- 
servative and be biased toward Incorrect exclusion of 
beats rather than Incorrect Inclusion of beats. One 
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suitable ECG classification methodology suitable for 
carrying out this invention is depicted in process 54 on 
FIG. 10. 

Next the SCG waveform 25 is cut up into wavelets 
5 based upon the position- in-time of the simultaneously 
recorded electrographic R waves. In general, the R wave 
peak is identified in the ECG 24 and a time window cor- 
responding to 100% of the measured R to R interval is 
applied to the SCG waveform. This window is asymmetri- 

10 cal around the R wave peak, cind data is collected 20% 
prior to the R wave peak and 80% after the R wave peak, 
as shown by wavelet windows wl through w6 in FIG. 6. 
Thus, the entire SCG waveform resulting from normally 
conducted beats is partitioned into wavelets. This 

15 process corresponds to 56 in FIG. 10. As shown in FIG. 
6 the process 54 excludes premature ventricular contrac- 
tion event (PVC) 7 and the corresponding SCG data are 
not folded into the wavelet series, shown in the drawing 
as 28-33. 

20 The wavelets etre next cross -correlated to ascertain 

their similarity. This corresponds to process 58 in FIG 
10. Self -similar wavelets are collected into sets. In 
general, most wavelets from a given data collection 
session fall into one class or set, called the dominant 

25 family- Classes with a small number of members (less 
than 25%) are deleted. This corxesponds to process 60 
on FIG. 10. 

The remaining sequence of retained (dominant) wave- 
lets are added cOid averaged together to form a canonical 

30 form wavelet which typifies the waveform taken from a 
specific patient during a particular phase of a par- 
ticular stress test. This canonical form wavelet is 
called the "representative waveform"^ for the particular 
case under consideration. This process is depicted 

35 schematically in FIG. 7 where wavelet type c dominates, 
with the largest membership. Line 34 represents the cut 
off for wavelet types . The sets below the line do not 
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have sufficient: membership to be considered useable sets 
so they are deleted. In the drawing type b, d and e 
wavelets are deleted, while wavelet types a and c are 
combined together to form the representative waveforms . 
5 At present the representative wavefojon of only one 

family is used for the advisory diagnosis . The wavelets 
wl through w6 are mapped into the histogram to show how 
the histogram is developed. At present any set with 
less than 25% of the total wavelets is excluded from the 
10 representative waveform. These steps correspond to 
processes 62 and 64 in FIG. 10. 

Although cross-correlation of data is well known it 
must be applied with care to the present data. 
Relatively brief corruptions of the data due to coughing 
15 result in a large niunber of families with very low 

membership, while veary heavy repetitive respiration can 
generate a smaller number of more heavily populated 
families. In general the correlation coefficients are 
varied to class and reclass the data until about 90% of 
20 the wavelets fall into a small (6) number of families. 
Data which doesn't fit this criteria is suspect, 
although the apparatus will still generate an advisory 
diagnosis . 
SCO Nomenclature 
25 FIG. 2 shows a schematic representation of simul- 

taneously recorded ECG and SCG data. The electrocar- 
diographic events depicted on the drawing correspond to 
the P wave 36, the QRS complex 37,38,39, and the T wave 
40. In general, the P wave corresponds to the contrac- 
30 tion phase of the atria of the heart which forces blood 
into the ventricles of the heart. The QRS complex 
arises from depolarization potentials generated by the 
contraction of the ventricular muscles of the heart. 
This contraction phase is followed by the T wave which 
35 generally represents the repolarization of the muscle 
tissue of the ventricle. 
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Through the use of simultaneously recorded ECG's, 
and echocardiogram data, certain repetitive points can 
be consistently identified on the SCG wave form set 
forth in the figure. For example The MC event 42 cor- 
5 responds to mitral valve closure. The AO event 43 cor- 
responds to the opening of the aortic valve. The RE 
event 44 corresponds to the rapid ejection of blood 
through the aorta. The AC event 46 corresponds to the 
aortic valve closure. The MO event 47 corresponds to 

10 the opening of the mitral valve. The RF phase 48 cor- 
responds to rapid ventricular filling. The AS 49 
inflection corresponds to the atrial systole. 
Identification of these wave form features on the fly is 
done based upon slope and amplitude information. The 

15 preferred decision rules applied to the representative 
waveforms are set forth as follows: 

The MC event is the first peak occurring in time 
after the peak of the electrographic R-wave. 

The AO event is the first valley after the MC. 

20 The RE event is the neact peaik after AO. 

The AC event is first peak after the end of the 
electrographic T-wave. T-wave are notoriously hard to 
locate due to the low slope. In the absence of success- 
ful T-wave detection based upon slope information the 

25 approximate T-wave position is defined based on measured 
R-R interval, the corresponding AC location is defined 
as a result. 

The MO event point is taken as the second valley 
following the AC event. 
30 The RF event is the next peak after MO. 

The AS event is taken as the last peak before the 
onset of the Q-wave in the ECG. 

The SCG waveform is scaimed according to these rules 
to locate these points or events on the SCG, and cor- 
35 responds to process 66 on FIG. 10. 
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The vaxious 1:line lnt:ervals between events as well 
their amplitudes and slopes are also measured, these are 
referred to as measured parameters. 

A number of quantitative parameters may be computed 
5 from these seismocsLrdiographic and electrocardiographic 
events. These parameters include: the Q to AC time 
interval, the Q to RE time interval. The RE to AC time 
interval^ the AC to RF time interval, the MC-AO ampli- 
tude, the AO-RE amplitude, the RF amplitude, the AS 
10 amplitude, as well as the MC-AO slope, the AO-RE slope 
and the MO-RF slope. These computed parameters are 
called the calculated parameters. 

These calculated parameters and measured parameters 
taken together are simply called the parameters. Certain 
15 of these quantitative parameters can be used to 

calculate traditional measures of cardiac performance, 
including the pre-ejection period (PEP), the left 
ventricular ejection time (LVET) , as well as the ratio 
of PEP to LVET. 
20 These and other non-traditional parameters can be 

used for developing the advisory diagnosis, at present 
the most useful parameter set known is set forth as the 
table of FIG. 8. The parameter set disclosed is all 
defined within substantially one beat of the heart. In 
25 other words, this illustrative and preferred parameter 
set is draw from underlying cardiac events occurring 
within one cycle of heart action. The formation of the 
parameter set corresponds to process 68 on FIG. 10. 
Neural Net Background 
30 A neural network is used to create the weighting 

coefficients used by the computation system 14 for the 
advisory diagnosis . 

Multi-layered neural networks are an emerging super- 
computer technology which utilizes leaming-by-back- 
35 propagation to categorize or classify certain patterns. 
This technology has been applied to bond rating, 
mortgage application evaluation, protein structure 
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determination, backgammon playing and hand-written digit 
determination. Neural network analysis can take the 
place of traditional statistical analysis of data to 
determine weighing coefficients, and is especially 
5 useful were training by example is preferred to 
programmed pattern analysis . 

AS used herein, a neural network is a system that 
produces an output vector that is a function of an input 
vector. The mapping fxinction between the input and 

10 output vectors is learned. In general, a neural network 
will include a number of processing elements, called 
neural units, arranged in layers. Interconnections are 
made between units in the successive layers. A typical 
network has an input layer, an output layer, and one or 

15 more hidden layers in between them. The hidden layer 
facilitates solutions to highly non-linear problems. 
Each of the processing units functions in a way which is 
analogous to a biological neuron. 

The original theoretical approaches toward under- 

20 standing neural networks were based upon the idea that 
when two neurons in the brain are active, the connection 
or synapse between them is strengthened. One early 
rule, developed by D. O. Hebb is described in his book, 
"The Organization of Behavior", Wiley, 1949. The 

25 Hebbian rule states that when two neurons are firing 
simultaneously, an association link between them is 
strengthened. Accordingly, the next time either of the 
two neurons fires, the other one is more likely to fire 
as well. Although the Hebbian rule is useful qualita- 

30 tively, it fails to eaqplain the learning process. For 
example, under a Hebbian rule, the coimection strengths 
between neurons grow without bound, or if maximums are 
placed upon connection strengths, these maximums are 
always reached. In 1962, Frank Rosenblatt introduced 

35 the Perceptron model in his book, "Principals of Neural 
Dynamics", Spartan, 1962. The Perceptron model includes 
input, hidden, and output layers consisting of 
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processing elements. An Input stimuli presented to the 
Input layer provides Information to the hidden layer* 
In a similar fashion, the hidden layer provides Informa- 
t tlon to the output layer* In the Perceptron models all 

5 the learning takes place at the output layer- Under the 
Perceptron model, many problems have been experimentally 
and mathematically shown to be representable by connec- 
tion strengpths between layers. 

In operation, the Perceptron model modifies the 
10 strengths of weighted connections between processing 
elements to learn an appropriate output vector or 
response, corresponding to a particular Input vector. 
The modification of the connection weights occurs when 
an Incorrect output response Is given. The modification 
15 of the weight Is done on an Interatlve basis, and 

changes the amount, or rate of transfer, of information 
from the Input layer to the output layer, so that even- 
tually, the appropriate output response will be 
provided. In retrospect, the Perceptron system is now 
20 understood as unable to properly adjust more than one 
layer of modifiable weights. Based upon this fact, it 
was predicted that no learning mechanism for a multiple 
layer system with modifiable weights was possible 
because none existed (Minsky & Pappert, 1969 in 
25 "Perceptron"). However, the more recent development of 
the "Back Propagation" algorithm as set forth in 
"Parallel Distributed Processing, the Mlcrostructure of 
Cognition", Rumelhart, Hinton and Williams, MIT Press, 
1986, has resulted in the development of multiple layer 
30 adaptive learning neural networks. Like the Perceptron 
model, the idea of errors is Introduced to correct or 
teach the machine to recognize a pattern. However, in a 
back propagation system, at each output processing 
element of the network, the error is quite easily 
35 realized, and this error is used to modify the strength 
of the connection between the processing element and the 
output processing element. This process permits the 
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alloca-tlon of an error value to hidden level elements, 
and permits the assignment of weight coupled to the 
hidden processing elements to be adjusted. However, 
most neural network algorithms, such as the back 
5 propagation algorithms, are computation intensive and 
require supercomputer capacity to efficiently explore 
pattern learning processes. However, once the weighing 
coefficients are determined, they may be utilized by 
substantially simpler computational machines to apply a 
10 learned response to new data sets or test cases. 
Computation and Use of the Weighing Coefficients 

The SCG data taken from a specific patient is called 
a "case". Fifty seven cases of patients with disease, 
and fifty seven cases of patients without disease were 

15 used to train a neural network. 

Each one of the fifty seven cases was prepared for 
the neural network by parameterizing each of the three 
representative waveforms collected for the case. The 
specific parameters formed are set forth in the table of 

20 FIG. 8. Thus each case was represented by 48 

parameters. This large group of parameters was supplied 
to a imiltilayered neural network simulator called 
BigMet, running a supervised learning back propagation 
algorithm on a Cray 2 supercomputer. A correction 

25 coefficient of (e=.l) resulted in mindLmal inter-itera- 
tion jumps in weighting. The output layer consisted of 
one node, producing a value of one, for the diseased 
state and a value of zero for the healthy state- As the 
result of training, the neural net selected weighing 

30 coefficients for the parameters as shown in the table of 
FIG. 9. 

The table of FIG. 9 reflects the relative importance 
of the various parameters. In general, coefficients 
which are large, correlate strongly with heazrt disease 
35 in the learning set. For example, parameters p5,p7 and 
pl4 are the most significant parameters of those tested. 

The optimal network for generating weighting coef fi- 
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clen-ts and the best parameters for presenting the cases 
are not intuitively obvious. Networks, for example, can 
over learn which results in 100% sensitivity and 
specificity on the learning set but very poor perfor- 
5 mance on the test set data. 

It is expected that weighting coefficients presented 
here will be modified and optimized as additional cases 
become available for analysis. It can also be 
anticipated that other effective parameters will be 
10 uncovered. For this reason the parameter choice and 

weighting values presented herein should not be taken as 
limiting but rather as illustrative of the invention • 

In use the apparatus collects the SCG and ECG data 
from the patient. Appropriate representative waveforms 
15 are formed for both the pre-exercise and the post exer- 
cise phases of the stress test and parameterized in 
accordance with PIG. 8. The parameters are multiplied 
by the appropriate weighting coefficient drawn from the 
table of PIG . 9 resulting in a "weighted parameter" . 
20 This corresponds to the process 70 on FIG. 10. 

The collection of the several weighted parameters is 
summed and multiplied by an determined "offset". The 
offset for the weighting coefficients presented in FIG. 
9 is ( w=.7241). The application of the offset cor- 
25 responds to the process 72 in FIG. 10. 

The summation and multiplication process results in 
the generation of a number which lies between unity (1) 
and zero ( 0 ) . Unity indicates disease with high con- 
fidence and zero represents disease free with high con- 
30 fidence. Intermediate values indicate varying degrees 
of indeteanninacy or "confidence". This process 
generates an "interpretive number", and corresponds to 
process 74 as shown on FIG. 10 • 

The theory and process of diagnostic stress test 
35 seismocardiography/ relies on comparisons between the 

pre-exercise or baseline SCG data with post exercise SCG 
data to determine if stress induced ischemia has altered 
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-the wall motion of the heart. Since the interpretation 
process is applied to pre-exercise and to post-exercise 
waveforms the underlying waveform comparison is 
inherently made. 
5 This interpretive number is mapped onto the 

diagnostic display as ''abnormal"(l-.7) , "borderline 
abnormal" {.?-• 6), "inconclusive'^t .6-,4) , "borderline 
normal "(.4-. 3) and "normal" ( .3-0.0) • This mapping opera- 
tion folds the confidence factor into the generation of 

10 the advisory diagnosis. This corresponds to the process 
74 on FIG. 10. Process 76 generates the hard copy shown 
in FIGS. 3, 4 and 5. The advisory diagnosis 78 is dis- 
played in one panel of the display |. as shown on FIG. 4. 
Twenty cases of SCG data taken from patients with 

15 coronary artery disease and twenty cases from patients 
without known coronary artery disease were used to test 
the initial neural network. Testing of the neural net- 
work is a process similar to the data interpretation 
processes carried out by the diagnostic system dis- 

20 closed. The testing process is carried out by para- 
meterizing cases that have not been used for training. 
The network applied the previously computed weighting 
coefficients to the parameterized data in a fashion 
analogous to the method performed by the diagnostic 

25 apparatus, and in a similar fashion the net an interpre- 
tive number. This advisory diagnosis is graded by com- 
pearing it with the known disease factors for the test 
cases. This grading is done on a "sensitivity " and 
"specificity" basis • 

30 Sensitivity is the ability to predict disease in 

patients with disease. It can be computed as the ratio 
of "true positive" to the sum of "true positive" and 
"false negatives". Specificity refers to the ability to 
predict wellness. It is taken as the ratio of "true 

35 negatives" to the sum of "true negatives" and "false 

positive". The best performances based upon neural net 
calculated weighting coefficients was 85% sensitivity 
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and 75% specl£lclt:y , which exceeds -the performance of 
-the best: human readers . 
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iTE CIiAXMs 

1. Seismocardiographic apparatus for collecting 

(16) cardiac information from a patient's body, and for 
5 displaying (18) an advisory diagnosis derived from said 
cardiac information to a diagnostician comprising: 

compression wave transducer means (10)/ adapted 
for coupling to the external surface of said patient's 
body (12)/ for collecting time series compression wave 
10 information resulting from the motion of said patient's 
heart/ and for producing seismic waveform data; 

segmenting means (56) for partitioning said 
seismic waveform data into a plurality of wavelets; 

cross -correlation means (58) for collecting 
15 similar ones of said wavelets into sets where set 

members have similar shape, and for generating a set of 
representative waveforms; 

event extraction means (66) for identifying 
event locations within said representative waveforms; 
20 parameter measuring means for collecting para- 

meter set data from said representative waveforms and 
for generating a peurameter data set (68); 

weighting means for applying a weighting 
coefficient to each corresponding psirameter of said 
25 parameter set and for generating a set of weighted para- 
meters (70) ; 

interpretation means for mapping said set of 
weighted parameters into an advisory diagnosis (72/ 74). 

30 2. A method/ performed by a machine, to form a 

preliminary and advisory diagnosis (76) of a disease 
condition comprising the steps of: 

a) collecting physiologic time series data 

(50); 

35 b) parameterizing said time series data which 

transforms said data into a set of parameters (68); 
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c) weighting said paraxnetiers by applying a 
predefined set: of weighting coefficients to each of said 
parameters, forming a set of weighted parameters (70); 

d) summing said weighted parameters to form an 
5 interpretive niimber indicative of the presence of 

disease underlying said physiologic time series data 
(72). 

3 • Apparatus for detecting a disease condition in 

10 a patient comprising: 

compression wave sensor means (10) for 
collecting a physiologic waveform representing time 
series of patient data; 

processing means (22) for parameterizing said 
15 time series data and for generating a seqfuence of para* 
meters ; 

weighting means (22) for applying a set of 
predefined weighting coefficients to said parameters and 
for generating a set of weighted parameters; 
20 computing means (22) for scaling and trans 

forming said weighted parameters into a diagnostic deci- 
sion space. 

4* A method of generating weighting coefficients 

25 (FIG. 9) for use in a seismocardiographic analysis 
system, comprising the steps of: 

collecting seismocardiographic data correspond- 
ing to a diseased state and corresponding to a non- 
diseased state; 
30 parameterizing said data, generating para- 

meterized data; 

presenting said parameterized data to a multi- 
layer neural network for back-propagation supervised 
learning; 

35 collecting weighting coefficients from said 

network after a fixed ntamber of network iterations . 
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5. A method of formulating an advisory diagnosis 

of a patient's disease state, comprising the steps of: 

a) collecting electrocardiographic and seis- 
mocardiographic time series data from said patient (SO); 
5 bl) psirameterizing said time series data by 

locating at a plurality of events, 

wherein said events are selected from the 
group consisting of: 

electrographic , R-wave events; 
10 seismocardiographic , MC events, AO events, 

RE events, PE events, AC events, MO events, RF 
events, AS events; 

b2) parameterizing said time series data by 
making a plurality of measurements, 

15 wherein said measurements are selected 

form the group consisting of: 

the Q to AC time, the PEP time, the LVET 
time, the ratio of PEP to LVET, the MC to AO 
amplitude, the MC to AO slope, the AO to RE 

20 amplitude, the AO to RE slope, the ratio of AO 

to RE slope divided by the MC to AO slope, the 
sum of the AO to RE slope and the MC to AO 
slope, the Q to RE time, the AC to Q time, the 
AC to RF time, the MO to RP amplitude, the MO 

25 to RT slope, the AS amplitude, 

c) weighting said parameters by applying a 
predefined set of weighting coefficients to each of said 
parameters, forming a set of weighted parameters (70); 

d) summing said weighted parameters to form 
30 an interpretive number indicative of the presence of 

disease underlying said physiologic time series data 
(72). 
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PARAMETERS 


MEASURED 
RESULTS 


NORMAL 
RANGE 


AORE / MCAO AMP 
AORE / MCAO SLP 


POSTEX/BASE = 0.617 
P0S7EX/BASE = 1.080 




PARAMETER ANALYSIS INDICATES ABNORMAL. 


BASEUNE 
SYSTOUC INTs 


MEASURED 
RESULTS 


MEASURED 
NORMAUZED 


NORMAL 
RANGE 


TOTAL SYS 

PEP 

LVET 

PEP/LVET 


460 MS 
116 MS 
344 MS 
0.34 


443 MS 
115 MS 
326 MS 


'409-439MS 
96-1 20MS 
304-330MS 
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SYSTOUC INTs 


MEASURED 
RESULTS 
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NORMAUZED 


NORMAL 
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300 MS 
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PI 


SYSTOUC INTFRVALi O TO AC 


1 1 19 


P2 


PRE— EJEiCTION PERIOD 


msec 


P3 


LEFT— VENTRICULAR EJECTIONI TIME 




P4 


PRE—EJECTION PERIOD / 
LEFT—VENTRICULAR FVJFCTIQN TIME 




P5 


MC TO AO AMPI minp 


9 


P6 


MC TO AO ^1 OPF 




P7 


AO TO RP AUPi m tnr 


9 


PR 


AH Tn PP CI ADP 




PQ 


An Tn PP CI noc* ✓ kir* Tn An oi nocr 




P10 


AO TO RE SLOPE + MC TO AO SLOPE 


g/msoc 


P11 


Q TO RE INTERVAL 


msec 


PI 2 


DIASTOUC INTERVAL: AC TO Q 


msec 


PI 3 


AC TO RF INTERVAL 
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PI 4 


MO TO RF AMPLITUDE 
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PI 5 


MO TO RF SLOPE 


g/msec 
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AS WAVE AMPUTUOE 
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6.9 


P6 


-1.4 


-8.8 


-3.6 


P7 


29.0 


-31.0 


-3.5 


P8 


-6.4 


15.4 


27.9 


P9 


-6.7 


4.6 


3.7 


PIO 


-5.4 


3.5 


16.5 


Pll 


-7.7 


-10.7 


6.0 
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7.9 


13.5 


-10.8 


P13 


-0.3 


-8.1 


1.4 


P14 


10.2 


19.2 


29.1 


P15 


-3.9 


-24.1 


-42.2 


P16 


5.4 


16.2 


-9.5 
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